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Globally Polarized Quark-gluon Plasma in Non-central A + A Collisions 
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Produced partons have large local relative orbital angular momentum along the direction opposite 
to the reaction plane in the early stage of non-central heavy-ion collisions. Parton scattering is 
shown to polarize quarks along the same direction due to spin-orbital coupling. Such global quark 
polarization will lead to many observable consequences, such as left-right asymmetry of hadron 
spectra, global transverse polarization of thermal photons, dileptons and hadrons. Hadrons from 
the decay of polarized resonances will have azimuthal asymmetry similar to the elliptic flow. Global 
hyperon polarization is studied within different hadronization scenarios and can be easily tested. 

PACS numbers: 25.75.-q, 13.88.+e, 12.38.Mh, 25.75.Nq 



Strong transverse polarization of hyperons has been 
observed in unpolarized p + p and p + A collisions since 
the 1970'sQ. Given the beam and hyperon momenta p 
and ph, hyperons produced in the beam fragmentation 
region are found transversely polarized in the direction 
perpendicular to the hyperon production plane, nn — 
P x Ph /\p x Ph I • Polarizations of A, H and H are negative 
while E and E's are positive. In the meantime, A and 
£1 are not transversely polarized. Although the origin 
for such striking transverse hyperon polarization is still 
in debate, one can relate it to the single-spin left-right 
asymmetries observed in hadron-hadron collisions with 
transversely polarized beam [2(, which in turn can be 
attributed to the orbital angular momenta (p.a.mj of 
the valence quarks in a polarized nucleon (HSU or 
fragmentation functions of transverse polarized quarks 
as well as twist-3 parton correlations in nucleons 0. 
It has also been suggested 0,0 that hyperon polarization 
could disappear due to the formation of QGP. 

In this Letter, we show that parton interaction in non- 
central heavy-ion collisions leads to a global quark polar- 
ization along the opposite direction of the reaction plane, 



n b =px b/\p x b\, 



(1) 



as determined by the nuclear impact parameter b. This 
global polarization is essentially a local manifestation of 
the global angular momentum of the colliding system 
through interaction of spin-orbital coupling in QCD. It 
will have far reaching consequences in non-central heavy- 
ion collisions, such as left-right asymmetry of hadron 
spectra in the reaction plane, global transverse polariza- 
tion of direct photons, dileptons and hadrons with spin. 
Within different hadronization scenarios, we will discuss 
hyperon polarization as a result of the global quark polar- 
ization. Possible contribution from final state hadronic 
interaction will also be discussed. 

Let us consider two colliding nuclei with the beam pro- 
jectile moving in the direction of the z axis, as illustrated 
in Fig.H We define the impact parameter b (along x) as 



the transverse distance of the projectile from the target 
nucleus and the reaction plane as given by fib (along y) 
in Eq. (I). Partons produced in the overlapped region 
of the collision will carry a global angular momentum 
along the direction opposite to the reaction plane (— y). 
A thermalized QGP requires final state parton interac- 
tion. Given the nature of partonic interaction at high 
energy, the global angular momentum would never lead 
to a collective rotation of the system. It will, however, be 
manifested in the finite transverse (along x) gradient of 
the average longitudinal momentum p z (x,y,b) per pro- 
duced parton. We assume for the moment that p z (x, y, b) 
is independent of the longitudinal position and is just an 
average value. Given the range of interaction Ax, two 
colliding partons will have relative longitudinal momen- 
tum Ap z — Axdp z /dx with o.a.m. L y ~ — AxAp z along 
the direction of fib. This relative o.a.m. L y will lead to 
global quark polarization due to spin-orbital coupling. 
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FIG. 1: Illustration of non-central heavy-ion collisions with 
impact parameter b. The global angular momentum of the 
produced matter is along — y, opposite to the reaction plane. 

The initial collective longitudinal momentum can be 
calculated as the total momentum difference between 
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participant projectile and target nucleons, whose trans- 
verse distributions (integrated over y) are, 



3.4 



Vn 



dx 2itRa 
+ [1 - (x zp fo/2) 2 ] arcsin 



r,2 

max 



(iTV2) 5 



1 - K 



2/n 



| ^/l- (2 + 6/2)2, ^i_ (i _6/2)2|( ? ) 



where i = x/Ra, b = b/RA and itU = 1.12A 1 / 3 is the 
nuclear radius in a hard-sphere distribution. 

Since the measured total multiplicity in A + A colli- 
sions is proportional to the number of participant nucle- 
ons ^3 j we can assume the same for the produced par- 
tons with a proportionality c(s). The average collective 
longitudinal momentum per parton is then, 



Pz(x,b) 



2c(s) 



dN^Jdx- 



dN^Jdx 



dN^Jdx + dN^Jdx' 



(3) 



The distribution p z {x, b) is an odd function in both x and 
b and therefore vanishes at x = or b = 0. As shown 
in Fig. [2] (as dashed lines) in units of Apo = Ayfs/2c{s) 
and as a function of x for different values of b, it is a 
monotonically increasing function of x until the edge of 
the overlapped region \x ± b/2\ — 1 where it drops to 
zero. The transverse gradient dp z /dx, shown as solid 
lines in unit of dpo/dx = \/s/2c(s)Ra, is an even func- 
tion of x and vanishes at b = 0. It increases almost 
linearly with b for small and intermediate values of b. 
Except the singular behavior at the boundary of the 
overlapped region which is caused by the assumed hard- 
sphere nuclear distribution, dp z /dx is approximately uni- 
form across the transverse cc-direction. The classical rel- 
ative o.a.m. is then L y = — (Ax) 2 dp z /dx for partons 
separated by Ace in the transverse direction. Due to 
transverse expansion, dp z /dx will decrease with time ac- 
cording to {R A - b/2)/[R A - b/2 + v x {b)(r - r )], where 
v x (b) is the transverse flow velocity in the reaction plane. 

In An + Au collisions at a/s = 200 GeV, the number of 
charged hadrons per participant nucleon is about 15 [13. 
Assuming the number of partons per (meson dominated) 
hadron is about 2, then c(s) ~ 45. Given Ra = 6.5 fm, 
dp Q /dx ~ 0.34 GeV/fm and L = -{Ax) 2 dp /dx ~ 1.7 
for Ace = 1 fm. 

We can relax the approximation of uniform distribu- 
tion of p z (x,b) in the longitudinal direction by iden- 
tifying pseudo-rapidity with the spatial rapidity rj — 
0.5 ln(t + z)/(t — z). According to experimental studies 
of hadron production in p + A and A + B collisions, the 
collective longitudinal momentum p z and dp z /dx are dis- 
tributed across a broad range of rapidity and peaks in the 
forward (backward) region for a given layer of dense mat- 
ter at positive (negative) x. The position of the peak in 
rapidity should increase with \x\. Note that even though 



p z (x,b) vanishes at around x — 0, dp z /dx and L y are 
still finite for b ^ as shown in Fig. [2J Averaging over 
the cc-direction will result in finite dp z /dx and L y around 
central rapidity region in non-central heavy- ion collisions. 





4 








3.5 


xT 




X 


3 


CL 






2.5 






~° 


2 


Cl 




"O 


1.5 


X 




\. 


1 








0.5 


cl 












0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x/R A 

FIG. 2: (dp z /dx)/(dpo/dx) (solid) and p z (x, fe)/4po (dashed) 
as functions of xj Ra for different values of b/ Ra- 

To study quark polarization due to parton collisions 
with a fixed direction of o.a.m., we consider quark scat- 
tering with fixed impact parameter ccV- For given initial 
relative momentum (E,p) and final spin A/2 of the quark 
along n;,, the cross section is, 



da 



A 



l\(q T ,k T ,E), 



d 2 x T ~" J (2tt) 2 (2tt) 5 
9 2 

Z\ = 2 ^ 2E y U X {p q )4{q T ){i> + m q )4(k T )u x {p k ), (4) 

within the screened potential model ^lj, where Ao(qx) = 
g/iOr + M 2 ) is the screened static potential with Debye 
screen mass Ct is the color factor associated with the 
target. Average over spin polarization is implied for the 
initial quark and p q (k) = P + fir (fa) is the final quark 
momentum. For small angle scattering, qT,kx ~ fJ> <C E, 



% » U Q (q T )A (k T ) 
9 2 



1 _ - A (gr - k T ) -(n b x p) 
1 2E(E + m q ) 



The first term gives the unpolarized cross section 

da da a- da_ „ „ , 

AC T a 2 s K 2 (iix T ), 
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(5) 



(6) 



or da/dq T = Ct^koq / (q 2 ^ + fi 2 ) 2 in momentum space. 
The second term gives rise to a polarized cross section 
dAa/d 2 XT = da + /d 2 XT — da^/d 2 XT, 

dAa p '■ (x T x ri b ) , _ 9 T ^ , . _ _ , . 



d 2 XT 



E{E + m q ) 



where cc<r = xt/xt and -fC n 's are modified Bessel func- 
tions. It is evident that parton scattering polarizes 
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quarks along the direction opposite to the parton reac- 
tion plane determined by the impact parameter xt, the 
same direction of the relative o.a.m. This is essentially 
the manifest of spin-orbital coupling in QCD. Ordinar- 
ily, the polarized cross section along a fixed direction nj, 
vanishes when averaged over all possible direction of the 
parton impact parameter xt- However, in non-central 
heavy-ion collisions the local relative o.a.m. L y provides 
a preferred average reaction plane for parton collisions. 
This will lead to a global quark polarization opposite to 
the reaction plane of nucleus- nucleus collisions. This con- 
clusion should not depend on our perturbative treatment 
of parton scattering as far as the effective interaction is 
mediated by the vector coupling in QCD. 

Averaging over the relative angle between parton xt 
and nuclear impact parameter b from —tt/2 to tt/2 and 
over xt, one can obtain the global quark polarization, 

P q = -TTfip/AEiE + mq) (8) 

via a single scattering for given E. If we take p = 
Ap z /2 — xt -xdpz/dx before averaging, the result will be 
similar but numerical evaluation is needed. In the limit 
m q = and p -C /i, one expects P q ~ —Ap z /fj,. Given an 
average range of interaction Ax -1 ~ fj, ~ 0.5 GeV and 
dpo/dx = 0.34 GeV/fm for semi-peripheral (b = Ra) 

collisions (see Fig. HJ at RHIC, P q 0.3. Multiple 

scattering will further increase the polarization. 

In nonrelativistic limit for massive quarks, m q 3> p, /j,, 

P q ^ -TTfip/8m 2 q . (9) 

In the same limit, the spin-orbital coupling energy is 
E LS = (L ■ S)(dV /dr)/(rm 2 q ). Given the range of 
interaction r ~ dVa/dr ~ — /i 2 and L ~ p/a 1 , 

Els I A* ~ "MP/ m q is j us t the above quark polarization. 

The global quark polarization opposite to the reac- 
tion plane will have many observable consequences in 
non-central heavy-ion collisions if an interacting QGP 
is formed. One expects to see left-right asymmetry in 
hadron spectra at large rapidity similar to the single- 
spin asymmetry in p+p collisions 0, Bl IS 0> • Thermal 
photons, dileptons and final hadrons with spin will be 
similarly polarized. Since hadrons from the strong decay 
of polarized resonances have angular distributions that 
prefer the direction perpendicular to the resonances' po- 
larization, they will result in an azimuthal asymmetry 
with respect to the reaction plane, similar to the asym- 
metry due to elliptic flow. In the following, we will dis- 
cuss global hyperon polarization since it can be easily 
measured through the weak decay [l^ . 

To demonstrate the robustness of the qualitative fea- 
tures of the predicted hyperon polarization due to global 
quark polarization, we consider several hadronization 
scenarios. We consider first hadronization via parton 



recombination. In this case, not only the spin of po- 
larized quarks but also the relative o.a.m. can con- 
tribute to the final hadrons' polarization. Given hadron 
size Ax, the classical estimate of the relative o.a.m. is 
L y ~ — (Ax) 2 dp z /dx. In the nonrelativistic quark model, 
however, constituent quarks in the ground state are all in 
the s-wave state. The contribution of o.a.m. to hadron 
polarization resides in the total angular momentum of the 
constituent quark. Over all, the effective polarization of 
a constituent quark should be proportional to that of the 
valence quark. Alternatively, polarization of constituent 
quarks can be similarly estimated, e.g. Eq. @, assuming 
massive constituent quarks as point-like particles. 

We can categorize recombination into exclusive qqq — > 
H and inclusive qqq — > H + X processes. The production 
cross section of polarized hyperons from the exclusive 
recombination can be written as 

°H = E I 2 ?3 3 I # A ) 1 2 92 93 Rql RqZ Rqs, (10) 

where R q = (1 + \P q )/2 is the quark polarization prob- 
ability. Extracting R H from a H = uhRh that only de- 
pends on quark polarization P q , one can calculate the re- 
combination probability Rh = Rr+Rh an< ^ * ne hyP eron 
polarization Ph = (R H — R\j)/Rh- Given baryons' 
SU(6) wavefunctions in the quark model 0], the polar- 
ization for strange (P s ) and non-strange quarks (P q ), one 
can obtain P A = P s , R A = 3(1 - P 2 ); 

= (4P, -P.- 3P S P 2 )/P S , P E = 3 - AP q P s + P 2 ; 
P~ = (4P S - P q - 3P 9 P 2 )/P H , R E = 3-AP q P s + P 2 : 
Pn - 2P S (5 + P 2 )/P f2 , P f2 -6(1 + P 2 ). 

If P s ~ P q in the most likely case, we have the same 
Ph = P q with R H = 3(1 - P 2 ) for A, E and 3. 

It is difficult to estimate hyperon polarization from in- 
clusive recombination. However, it could be the domi- 
nant process for the bulk hadron production, especially 
for large values of P q . It is also required by entropy con- 
servation. The polarization of produced hyperons would 
be smaller than in the exclusive recombination but should 
be proportional to P q . 

The extreme limit of inclusive recombination is frag- 
mentation of polarized quarks, q — * H + X. Lon- 
gitudinal polarization of hyperons in a similar process 
e + e~ — > Z° — > qq — > A + X has been measured 0] 
and can be explained [lij by assuming that polarized hy- 
perons contain the initial polarized leading quark in its 
SU(6) wavefunction. One can similarly calculate hyperon 
polarization from fragmentation of transversely polarized 
quarks and obtain, 

Pa = n s P s /{n s + 2/ s ), P s = {Af s P q - n s P s )/3{2f s + n s ), 
P~ = {4n s P s - f s P q )/3(2n s + f s ),P n =P s /3, 
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where n s and f s are the strange quark abundances rela- 
tive to up or down quarks in QGP and quark fragmen- 
tation, respectively. Including production of unpolarized 
hyperons in the process, the effective hyperon polariza- 
tion from fragmentation of polarized quarks should in- 
crease with the fractional momentum zjj = pu/Pq- 

Because of the complication from hadronization, we 
cannot provide a model independent estimate of the fi- 
nal hyperon polarization. If quark recombination is the 
dominant hadronization mechanism, hyperons' polariza- 
tion will be determined by the quarks' polarization be- 
fore hadronization. According to our estimate following 
Eq. ijSJ , this would be in the order of tens of a percent and 
the hyperons' polarization would be around the same or- 
der. We can also provide other qualitative predictions of 
the global hyperon polarization Pu in non-central heavy- 
ion collisions: (1) Hyperons and their anti-particles are 
similarly polarized along the same direction perpendicu- 
lar to the reaction plane in non-central heavy-ion colli- 
sions. (2) The global hyperon polarization Pjj vanishes 
in central collisions and increases almost linearly with b 
in semi-central collisions. (3) It should have a finite value 
at small px and in the central rapidity region. It should 
increase with rapidity and eventually decreases and van- 
ishes at large rapidities. (4) High px ^> Ap z ~ \iLq 
quarks should not be polarized by parton scattering in 
the medium. However, hyperon polarization should per- 
sist at moderate px where recombination of unpolarized 
high pt quarks with polarized thermal quarks may still 
dominate. The polarization can be estimated within a 
recombination model. (5) Since hyperon's production 
planes are randomly oriented with respect to the reac- 
tion plane of heavy-ion collisions, the observed hyperon 
polarization in p + A collisions should not contribute to 
the global polarization as we have defined here, except at 
large rapidity region where directed flow is observed 0] . 
In this region, the non-vanishing (p x ) can provide an av- 
erage production plane nu for hyperons. According to 
the observed polarization pattern in p the global 

Ph will be enhanced for A, S, and 5 and reduced for 
£ in large rapidity region. In addition, one also expects 
Pa > P^ due to the directed flow. 

In principle, hadronic interaction with a given direction 
of relative o.a.m. can also lead to global hyperon polar- 
ization. Neglecting hyperon production in the hadronic 
phase, the dominant hadronic processes involving hyper- 
ons will be hyperon-7r scattering. The scattering ampli- 
tude in general involves both scalar and vector channels 
|l6| . One can show that global hyperon polarization from 
the vector channel is along while polarizations from the 
scalar channel and all the interferences are against the 
global quark polarization. Therefore, polarizations from 
different channels in hyperon- 7r scattering partially can- 
cel each other. In addition, the transverse gradient of 
the longitudinal flow should be significantly reduced in 
the hadronic phase due to prior strong transverse expan- 



sion in the reaction plane as demonstrated by the large 
elliptic flow measured at RHIC. Therefore, the net ef- 
fect of hadronic interaction should be small and should 
not change the final hyperon polarization significantly, 
in particular at large pr- We will leave detailed study to 
future work. 

In summary, produced partons are shown to have large 
local relative o.a.m. in non-central heavy-ion collisions if 
quark-gluon plasma is formed. Parton scattering with 
given relative o.a.m. is shown to polarize quarks along 
the same direction due to spin-orbital coupling. Such 
global quark polarization has many measurable conse- 
quences in high-energy heavy- ion collisions. Within dif- 
ferent hadronization scenarios, we predict that hyperons 
will be polarized along the opposite direction of the reac- 
tion plane. Effects of hadronic interaction are expected 
to be small and would not change the qualitative feature 
of our prediction. 
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Erratum 

After the publication of the manuscript, we found a 
factor of 2 missing from Eqs. (8) and (9). Therefore the 
global quark polarization in Eq. (8) should be 

P q = -khp/2E{E + m q ). (11) 

In the non-relativistic limit for massive quarks, m q <C 
p,fj,, the polarization becomes [Eq. (9)], 

P q « -7Tfip/4m 2 q . (12) 



We also made an error in the numerical estimate of the 
average longitudial momentum difference within a typical 
interaction range in non-central Au + Au collisions 
at RHIC and consequently the quark polarization follow- 
ing Eq. (8). The discussion following Eq. (8) should be 
replaced by: 



In the limit m q = and p 3> /x, one expects P q = 
— 7r/z/2Ap z . Given dpo/dx — 0.34 GeV/fm for semi- 
peripheral (b = Ra) collisions at RHIC and an average 
range of interaction Ax -1 ~ fj, ~ 0.5 GeV, Ap z ~ 0.1 
GeV, which is much smaller than the average transverse 
momentum transfer fi. Therefore for a reliable estimate 
of the quark polarization at RHIC, one needs to go be- 
yond the approximation of small angle scattering. 



